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We study the diffusion equation with a position-dependent, power-law diffusion coefficient. The 
equation possesses the Riesz-Weyl fractional operator and includes a memory kernel. It is solved 
in the diffusion limit of small wave numbers. Two kernels are considered in detail: the exponential 
kernel, for which the problem resolves itself to the telegrapher's equation, and the power-law one. 
The resulting distributions have the form of the Levy process for any kernel. The renormalized 
fractional moment is introduced to compare different cases with respect to the diffusion properties 
of the system. 
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I. INTRODUCTION 



Diffusion processes are usually described in terms of either differential or fractional equations which contain a 
constant diffusion coefficient. In many physical problems, however, that coefficient depends on the position variable 
and such dependence is important [Ij]. As a typical example can serve the transport in porous, inhomogeneous 
media and in plasmas. Modelling the aggregation of interacting particles must take into account nonlocal effects, 
since the particle mobility depends on the average density the coalescence of particles results from long-range 
interactions (the Poisson-Smoluchowski paradigm) and the corresponding evolution equations contain a position- 
dependent coefficient. That modelling can be accomplished directly, via the non-local Fokker-Planck equation, in 
which the term with the space derivative is multiplied by a kernel and integrated over the position A similar 
method, applicable to the Levy processes, consists in the integrating over the Levy index, with some kernel (the 
distributed order space fractional equation) [4]. The spatial inhomogeneity can be also taken into account as an 
external potential which may substantially change the diffusive properties of the stochastic system, in particular of 
the Levy flights 0. 

The Levy distributions constitute the most general class of stable processes and the Gaussian distribution is their 
special case. One can expect that the Levy (and non-Gaussian) distributions emerge in transport processes for which 
the observable values experience long jumps, e.g. due to the existence of long range correlations. The theory of 
the Levy flights is applicable to problems from various branches of science and technology. Moreover, the handling 
of specific and realistic systems often requires taking into account both memory effects and inhomogeneity of the 
media. As a typical example of the nonhomogeneous problem can serve the diffusion in the porous media; they often 
display the fractal structure and the diffusion on the macro- and mesoscale can be expressed by a stochastic equation 
driven by the Levy process In general, the transport in fractal media can be described by the fractional Fokker- 
Planck equation with a variable, position dependent, diffusion coefficient 0, The Levy flights bring about the 
accelerated diffusion in the reaction-diffusion systems ^IflJ and the probability distribution for that process is expressed 
by the fractional Fisher-Kolmogorov equation. The Levy processes are typical for problems of high complexity, in 
particular in biological systems 11] where the fractal structures are also encountered. For example, the lipid diffusion 
in biomembranes has the characteristics of the Levy process but it can no longer be regarded as Markovian. The 
theory of Nonnenmacher takes into account the memory effects, as well as the fractal structure of the medium; 
the diffusion coefficient depends on the variable diameter of the holes in the solvent through which the molecules 
jump. The application of the Levy processes is natural also in many social and environmental problems. Recently, it 
has been established [ll| that the people mobility, estimated by the bank notes circulation and studied in terms of 
the stochastic fractional equations, strongly depends on geographical and sociological conditions. Therefore, its study 
requires including position-dependent quantities. That problem is directly related to the spread of infectious diseases. 
It has been demonstrated in the example of the SARS epidemic and by means of percolation model simulations (T^ 
that the disease can spread very rapidly. Usually one assumes that the infection probability at a given distance is 
Levy distributed due to the long-range interactions but the process is local in time [15|. On the other hand, the 
percolation model of the epidemics developed in Ref.fl6l] is restricted to short-range interactions (is local in space) 
but it introduces the incubation times which obey the Levy statistics and then the model is non-Markovian. 

In Refs.firl. [Tsj . the master equation for a jumping process, stationary and Markovian, has been studied. That 
process is a version of the coupled continuous time random walk (CTRW), defined in terms of two probability 
distributions: the Poissonian waiting time distribution with the position-dependent jumping frequency and a jump- 
size distribution. The standard technique to handle such master equations is the Kramers-Moyal expansion which 
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produces the Fokker-Planck equation for the Gaussian jumping size distribution and it yields correct results for large 
times and large distances [19|. For the Levy distributed jumping size, the Fokker-Planck equation becomes the 
fractional diffusion equation, with the Riesz-Weyl fractional operator and the variable coefhcient D(x). Formally, it 
can be derived from the master equation by taking the Fourier transform and by neglecting the higher terms in the 
wave number expansion of the jumping-size distribution (the diffusion approximation) ^20] . The equation reads 

dt d\x\t^ ' ^ ' 

where 1 < /i < 2. Since the diffusion coefficient is just the jumping frequency, the medium inhomogeneity enters the 
problem via the a;-dependent waiting time distribution. For the Gaussian case (/i = 2), all kinds of diffusion, both 
normal and anomalous, are predicted [29| and they are determined by the jumping frequency. 

The Eq.([T]) can be regarded as a special case of a more general problem than the random walk and which traces 
back to the microscopic foundations of the nonequilibrium statistical mechanics. The well-known achievement of 
Zwanzig [ 2l| was the derivation of the non-Markovian kinetic equation for the probability distribution in the space 
of macroscopic state variables. More precisely, by starting from the Liouville-von Neumann equation for the density 
operator p: ihdp/dt = [H,p\, where H is the Hamiltonian of the system, one can obtain the generalized master 
equation: 

- dt'4>{t t') Y^WMt') - FMt% (2) 

where denotes the diagonal elements of the density matrix and F^j are the transition rates ^231 . Then the equation 
is non-Markovian and it contains the memory kernel 4>{t). Markovian equations like ([T]) follow from the generalized 
master equation if the memory effects are negligible. However, usually this is not the case. We have already discussed 
the examples of the Levy processes, with power-law tails of the distribution, which exhibit memory effects. In fact, 
the importance of these effects was realized a long time ago, e.g. in the description of the resonance transfer of the 
excitation energy between molecules [l^] . The detailed calculation for the anthracene molecules shows that the memory 
kernel is exponential and the generalized, nonlocal in time, master equation must be applied to get proper results 
for small times. One can expect that memory effects are still more pronounced for systems with the characteristic 
decay rate slower than exponential, that often happens for atomic and molecular systems [1^. Stochastic dynamical 
processes are generally nonlocal in time due to the finite time of the interaction with the environment. Moreover, for 
a stochastic system which is coupled to a fractal heat bath via a random matrix interaction [24] , the finite correlations 
emerge and its relaxation has to be described in terms of the generalized, non-Markovian Langevin equation, with 
the memory friction term [2^, [2^ . Also the speed of transport is affected by the memory. In the non-Markovian 
CTRW processes it hampers the dynamics and such systems are subdiffusive [131. Such processes are described by 
the generalized master equation, with a memory kernel, if the waiting time distribution possesses long, algebraic tails. 
That equation follows directly from the generalized Chapman-Kolmogorov equation which determines the probability 
distribution in the phase space [H, [2^ . 

By applying the nearest-neighbour approximation on the transition rates and taking the continuum limit [2^ , 
one obtains from the Eq.® the non-Markovian Fokker-Planck equation. In the presence of long-range correlations, 
however, the nearest-neighbour approximation is no longer valid. If the transition rates are symmetric and distributed 
according to the Levy statistics in the continuum limit, the Kramers- Moyal expansion produces the fractional deriva- 
tive, instead of the Gaussian. Then, for the variable diffusion coefficient D{x), the equation which corresponds to the 
Eq.([5]) becomes 



dp^{x, t) 
di 

where the operator 



/ K^{t~t')C^[p^{xX)]dt\ 
Jo 



(3) 



^^-K^^nix) (4) 

acts only on the x variable. The parameter 7 measures the rate of the memory loss. 

The Eq.Q is of interest both from quantum and classical point of view. In the atomic and molecular physics e.g. 
in a few- modes spin boson model [sot and the random-matrix theory [3l|, where the decay rate is slow, an equation 
analogous to the Eq.Q can be applied. The operator is then expressed in terms of the " superoperator" which 
represents an instantaneous intervention of the environment over the system (23j and it can assume a quite general 
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form. In the classical context, the Eq.Q has been discussed in Ref.[32']; the operator has the Fokker- Planck form 
in this case, with the constant diffusion coefficient and a potential force. 

In this paper we study the diffusion problem for systems which are driven by the Levy distributed transition rate 
and for which both the medium inhomogeneity and the memory effects are important. We assume D{x) — \x\~^ 
[9 > —1). The power-law form of the diffusion coefficient has been used to describe some physical phenomena, e.g. 
the transport of fast electrons in a hot plasma [s^ and the turbulent two-particle diffusion (s^l- It is also used in 
theoretical analyses of the fractional Fokker-Planck equation [sl, [s^, [s^, Hg], e.g. as an ansatz for the problem of 
diffusion in the fractal media 7, 8, 9, 39]. Obviously, for the Markovian case K^{t) = 5{t), the Eq.Q resolves itself 
to the Eq.Q. 

In Sec. II we solve the fractional telegrapher's equation which follows from the Eq.® for the case of the exponential 
memory kernel K^{t). The solution for an arbitrary kernel, expressed in the form of the Laplace transform, is derived 
in Sec. III. Moreover, the case of the power-law kernel is solved in details. In Sec. IV we derive the fractional moments 
and discuss their application to the description of the diffusion process. The results presented in the paper are 
summarized in Sec.V. 



II. THE EXPONENTIAL KERNEL 



If the memory effects are weak, we can assume that the kernel K^(t) decays exponentially. Then let us consider 
the following kernel: 

K^{t) = 76-^* (7 > 0) (5) 

which becomes the delta function in the limit 7 — > oo (the Markovian case). In this case, the integral equation, Eq.Q, 
reduces itself to a differential equation. It can be derived by inserting the Eq.® to the Eq.Q and by differentiating 
twice over time, in order to get rid of the integral. Finally, we get the following equation 

d^p^{x,t) dp^{x,t) d''[\x\-^pj{x,t)] 

dt^ +^^^=^^ — dw — ' 

which is a generalized and fractional version of the well-known telegrapher's equation. Originaly, the telegrapher's 
equation, which is the hyperbolic one, has been introduced into the theory of the stochastic processes by Cattaneo 
[40| in order to avoid infinitely fast propa gati on for very small times. Its fractional generalization describes e.g. a 
two-state process with the correlated noise Wr\ and it predicts the inhanced diffusion in the limit of long time. On the 
other hand, in the case of the divergent second moment, the telegrapher's equation with the Riesz-Weyl derivative 
results from the fractional Klein-Kramers equation for the Levy distributed jumping size [291] . In that equation, the 
parameter 7 has a sense of the damping constant in the corresponding Langevin equation. 

In the diffusion limit of small wave numbers, the Markovian equation, Eq.([T]), is satisfied by the Fox function 2 
[20I ]. Since our main objective is to study the diffusion problem, we restrict also the present analysis to that limit. 
We will try to find the solution of the Eq.([6|) in the same form as for the Markovian equation. Therefore we assume: 



p^{x,t) = NaH2'2 



a\x\ 



{ai,Ai), (a2,^2) 
(61, Bi), (62,^2) 



(7) 



where the time dependence is restricted to the function a{t) and N is the normalization constant. The method of 
solution, described in Ref . [l^] , consists in the inserting of the Fourier transform of the expression ([7]) into the Fourier 
transformed Eq. © . Then we determine the coefficients of the Fox function by demanding that the Eq. ^ should be 
satisfied in the diffusion limit, i.e. for small wave numbers. In fact, the latest assumption does not introduce any 
additional idealization since the equation itself is valid only in the diffusion limit. 
We start with the Fourier transform of the Eq.®; it reads 

d^p^ -f j^Mii^ = -K'^j\km\-\-%{x, t)]. (8) 

The Fourier transform of the Fox function can be expressed also in terms of the Fox function (for the definition 
and some useful properties of the Fox functions see Ref. [13] and references therein). Due to the multiplication rule, 
the product \x\~^p^{x,t) is the Fox function as well. Both sides of the Eq.® can now be expanded in series of 
fractional powers of \k\. We can satisfy the Eq.® by a suitable choice of the parameters of the function ([7]) and by 
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neglecting the terms higher than 
rhs, respectively, reads: •p^(k,t) « 1 - 



The results are the following. The expansion of the functions on the Ihs and 
iV/i^a~^|fc|^ and jr[|a;|~*p-y(a;, t)] « Nh^^ ^a^ ^ with the following coefficients: 



h^^^ = + e)/{2 + 9) and h,, = 
the order less than fj, is the necessary condition to satisfy the Eq. 



r(— /i)r(/i + 6) cos(/i7r/2) sin(^i|7r). The vanishing of all other terms of 



^ 2+9 

The solution takes the form 



Pj{x, t) = A^ai?2,' 



1,1 



a\x\ 



(I - 1^ -J-] (1-1^ 



^' 1)' (1 2+e ' 2+e) 



(9) 



and the coefficients bi and Bi are responsible for the distribution behaviour near x — 0. bi and Bi cannot be 
determined in the diffusive limit and they are meaningless from the point of view of the diffusion process; the values 
9 and 1 we have assumed correspond to the jumping process, considered in Ref . [20| . Generally, the Eq.Q is satisfied 
by the function ^ for any choice of the coefficients bi and i3i > 0, such that 6i ^ and i?i ^ 1 for ^ 0. The 
normalization factor N can be determined in a simple way from the formula N = [2x(— where xi~s) is the 
Mellin transform of the Fox function. A simple algebra yields 



N : 



r(i + 6i)r 



Alternatively, since |fc| is small, we can express the Fourier transform of the solution as 

P-y{k,t) « l-a^lkl" « exp(-(T^|/c|^), 



where 



(/i + i 



2r(-/i)r(/i + 0) cos W2) sin(f±f ) 



(10) 



(11) 



(12) 



The Eq. ((TT|l means that the solution of Eq.Q coincides with the Levy process in the limit k 0. Then the solution 
([9|) can be expressed in the form which is generic for any symmetric Levy distribution (4^ : 



(l-l/M,l//i),(l/2,l/2) 
(0,1), (1/2, 1/2) 



(13) 



The formula establishes the relation between the solutions ^ and Those expressions are equivalent only 

in the limit fc ^ and they behave differently for small \x\. We will demonstrate in Sec. Ill that the Levy process is 
the solution of the Eq.Q for any kernel. Therefore, the form (|13p is quite universal and we apply it in the following. 
The problem is reduced in this way to evaluating the fimction (j(t). 
Now, the Eq.® becomes the ordinary differential equation: 



7 " 



(14) 



where £^{t) = a~'^. We assume the following initial conditions: ^(0) = ^(0) = which correspond to the condition 
p^{x,0) = S{x). The Eg. pi)) has the structure of the equation of motion with a "friction term", a positive "driving 
force", and a "mass" I/7. The meaning of the quantity ^, the time evolution of which the Eq.([T4]) describes, remains 
to be determined. The variable ^, as well as ^, rises with time and finally the balance of "forces" , given by the equation 



is reached. Note that the above expression is equivalent to the Eq. (fT4jl in the Markovian limit 7 
p-yix, t) — Poo{x, t) — p{x, t) for t — > 00. The solution of the Eq. (fT5|) produces the result 



(15) 



00. Therefore 



a{t) = 



h. 



-i -i/(M+e) 



{t 00) 



(16) 
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FIG. 1: (Color online) Probability distributions for the case of the exponential kernel with fi = 1.5, calculated from the Eas. (|14p . 
(fT2]). and (O, for t = 50. The initial condition is Pj{x,0) = 5{x). 



which corresponds to the exact solution (for arbitrary time) for the Markovian hmit, p{x,t). 

The case of the constant diffusion coefficient, = 0, is a particular case and it can be solved easily. The solution of 
the Eq. pl]) leads to the following result 



ait) - 1 
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(17) 



For 9^0 and arbitrary time, the Eq. (jl4p can be solved by the numerical integration and the distribution p-yix, t) 
determined from the Ea. (jl3p . To evaluate the Fox function we use the general formula for its series expansion and 
then the Eq. (fT3| can be expressed in the computable form: 



^ LXJ 



r[l + (2n+l)/M] 
(2n+ 1)!! 



(f ) 



2n 



(18) 



Fig.l presents some exemplary probability distributions, so chosen to illustrate the influence of the memory on the 
time evolution. Since the serie s (|18p is poorly convergent, the evaluating of the distribution for large |x| required the 
quadruple computer precision |43j . The case with 7 = 3 is close to the Markovian one; a comparison with the case 
characterized by the long memory shows that the spread of the distribution slows down with the decreasing value 
of 7, i.e. for stronger memory (larger "inertia" in the Ea. ((T4)) ). In the limit t ^ 00 the curves which correspond to 
different 7 values and the same - coincide. 



III. THE GENERAL CASE 



The description by means of the Eq. ^ with the exponential memory kernel does not apply to systems with long-time 
correlations and small decay rate. In the study of realistic systems one encounters a variety of forms of the kernel; some 
of them are very complex. It is typical for natural signals that they do not represent a simple kinetics, characterised 
by a unique Hurst exponent. Random processes which take into account the whole spectrum of the time-dependent 
Hurst exponents serve then as useful models. This concept, applied to the fractional equations formalism, leads to the 
integration over the order of differentiation (the distributed-ordcr diffusion equation) A] and the kernel assumes the 
integral form: J f{a)t~"da. Reactions in polymer systems are also described by using complicated kernels (ij. (isj. 
Therefore, in this Section we consider the Eq.([3]) for the case as general as possible. We will demonstrate that the 
solution in a closed form can be obtained for the arbitrary kernel. 

The Eq.([3]) has the structure of the Volterra integro-differential equation with the kernel which depends on the 
difference of its arguments. Therefore, methods using Laplace transforms are applicable. Following Sokolov [32] . we 
apply a method of the integral decomposition which allows us to express the required solution by the solution of the 
corresponding Markovian equation ([T]). Let us define the function T(t, i) by its Laplace transform: 

T^(r,.)^^exp(-.^). (19) 
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If we know the function T(r, t), the probabihty distribution p^(x, t) can be obtained by a simple integration: 

POO 

p^{x,t) = / p{x,T)T{T,t)dT. (20) 



However, since the inversion of the Eq. (|T9)l is difhcuh for any kernel, it is expedient to get rid of T. To achieve that, 
we take the Fourier transform from the Eq. (|20p and eliminate the function T{T,t) by using the definition p9|l . The 
final solution is of the form of the following Fourier-Laplace transform 

p*Jk,u)^^p*(k,^\. (21) 



The above formalism can be applied for any kernel and any operator Cx] the main difficulty consists in inversion 
of the Laplace transforms. 

First of all, we find that if the Markovian process p{x, t) is Levy distributed in the diffusion limit fc — > 0, then the 
non-Markovian process is also the Levy process in this limit. Indeed, the Fourier transform of the Markovian solution 
is given by the Eq.([TT|), where a{t) follows from the Eqs. p^ and Then we take the Laplace transform from 

that expression and insert the result to the Eq. (PT|) : i5*(fc, u) = 1/u — F* {u)\k\'-'- . Finally, the inversion of the Laplace 
transform yields 

p^ik,t) = l-F{t)\kr, (22) 

which is just the Fourier representation of the Levy distribution for small To get the function F{t) we need to 
invert the Laplace transform 

^*(") = ^[«-TW) (23) 

and we assume that this inverse transform exists. The solution is given by the Eg. p^ . where a{t) — [F(t)]^^f^. 

We will consider two particular cases in detail. In the case of the exponential memory kernel ([5]) , discussed already 
in Sec. II, we have K*{u) — '-f/{u + 7) and the Eq. (PT|l produces the following result 

p;{k,u) = i -aor(l + «)7"|fcr"~^"+'Hw + 7)"", (24) 
' u 

where = Nh^^°'{K^h'^Q^ a/ ^)°' and a = /j,/(/i + 6). The above expression cannot be inverted in closed form. 
However, if we are interested only in large times, the last term in the Eq. ([M)) can be expanded around u — Q: 
{u + 7)^" ~ 7^" — a7^^"^^^u. Then the inversion of the Laplace transform yields 

p^(fc,i) = 1 -ao|fc|^ - ^t""^^ {t^oo) (25) 

and this expression demonstrates how the solution p^ approaches its asymptotic, Markovian form. The final solution, 
valid for large t, is given by the Eq.dni) with a = [ao(i" - ■yi""^)]^/''. 

The other physically important kernel has the power-law form, with long tails: 

K^{t) - f(T^ (0 < 7 < 1). (26) 

The equation Q with the kernel (j26|l is usually presented as the fractional equation with the Riemann-Liouville 
derivative p6| - which is equivalent to the Caputo operator for a special choice of the initial conditions - in the 
following form 

dpj{x,t) _ 1 



dt 



= oDr i:-APfix,t)]. (27) 



The power-law kernels are used to describe subdiffusive relaxation e.g. in the framework of the CTRW [23|. They 
emerge also as a result of the coupling to the fractal heat bath via the random matrix interaction [23|. To solve 
the Eq.Q we follow the same procedure as for the exponential kernel. The Laplace transform of the Eq.ipS)) reads 
K*{u) = u''-^ and the Eq.jn]) takes the form 

p*Jk, u)^-~ aor(l + a)7"|fcru-2-+7"-i. (28) 
' u 
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The inversion can be easily performed: 

^^^^'^^ = ^ r72r-7aV"l) '^'''^'°''" ^ ^ " ^^^^ 

Clearly, the above solution does not converge with time to the Markovian asymptotics, F{t) ^ t", in contrast to the 
case of the exponential kernel. 

To conclude this Section, we want to mention yet another approach to the Eq.Q, which is a direct generalization 
of the method applied for the telegrapher's equation in Sec. II. Inserting the expansion of the functions p^{k,t) and 
!F\\x\~^p^{x^ t)] to the Eq.Q confirms the finding that the solution can be expressed in terms of the Fox function H^'l^ 



and it is Levy distributed. The resulting equation is a generalization of the Eq.([T4]) and it determines the function 

^ = i^M:^ f KJt- t')r"''dt'. (30) 
dt Jq 

Mathematically, the Eq.((30| has the form of the nonlinear Volterra integro-differential equation. Since the numerical 
inversion of the Laplace transforms is not always an easy task (methods are often unstable), the numerical solving of 
the Ea. ((5Dl) could be a useful alternative to the Eq. 



IV. DIFFUSION 



The diffusion process is usually characterized by the time dependence of the second moment of the probability 
distribution: if this dependence is linear in the limit of long time, the diffusion is called normal. There are many 
examples of physical systems for which the variance rises faster than linearly with time (hyperdiffusion), or slower 
(subdiffusion) . Such behaviours are typical for transport in the disordered media [13] and systems with traps and 
barriers. In the realm of dynamical systems, a substantial acceleration of the diffusion is caused by the presence of 
regular structures in the phase space [48| . On the other hand, the subdiffusion appears in the non-Markovian version 
of CTRW, as a result of a non-Poissonian, power-law form of the waiting time distribution [23|. 

When we enter the field of the Levy processes, the situation becomes more complicated. The stochastic variable 
performs very long jumps and their size is limited only by the size of whole system. As a result, the second moment, 
as well as all moments of the order /i or higher, is divergent. Mathematically, that follows from the fact that the tail of 
the Levy distribution is the power-law: ~ Therefore, one cannot describe the diffusion process in terms of 

the position variance and some other quantity which could serve as an estimation of the speed of transport is needed. 
One possibility is to consider still the second moment but with the integration limits restricted to a time-dependent 
interval (the walker in the imaginary growing box) [49|. On the other hand, one can derive the fractional moments of 
the order d < ^. 

By the derivation of the moments of the probability distribution p^(x,t), Eg. ljlSp . we utilize simple properties of 
the Mellin transform from the Fox function: 

{\x\^) = 2 / x^p.y{x, t)dx = -a^x{-6 - 1) = -(T^y{&)y{l - -) sin(<57r/2). (31) 

Let us consider two quantities: the renormalized moment of order /z, defined by the following expression 

TW^ = lim e(|a;|'^-^) = -a'^r(//)sin(/i7r/2), (32) 

where we applied the property: r(a;) —fl/x for x — > 0, and then the fractional diffusion coefficient V''^\t) = 
In the Markovian case, defined by the Eq.((T]), the coefficient 2?*^^^ is useful to classify the diffusion: for 

6* < it rises with time, for 6* > it falls, and it converges to a constant for 6 — Q 20]. That pattern is consistent 
with the diffusion properties, defined in the ordinary sense, of the Fokker-Planck equation (/i = 2). Therefore, in the 
following wc will name all kinds of the diffusion - the subdiffusion, the normal diffusion, and the superdiffusion - 
according to the time dependence of the coefficient P*^^^ . 

We begin with the case of the exponential kernel. First we realize that, since = Nh^£,, the renormalized moment 
M'^ is directly related to the variable ^: = ;|A^/i^r(/i) sin(^7r/2)^. Therefore, the interpretation of the Eq. fH]) 
is straightforward: it describes the deterministic time evolution of the moment M.*^ . The diffusion properties of the 
system remain unchanged, compared to the Markovian case, because in the limit i — > oo both solutions coincide. 
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FIG. 2: (Color online) The fractional diffusion coefficient for the case of the exponential kernel with 7 = 0.05 (solid lines) and 
7 = 3 (dashed lines), as a function of time, obtained from numerical solving of the Ea. (ll4p . Results for three values of 9 are 
presented: 9 = —0.2 (upper lines for large t), 6 = (middle lines), and 9 — 0.4 (lower lines); fj, — 1.5. 

However, at small time the influence of the memory, which hampers both the spread of the distribution and the 
relaxation to the Markovian asymptotics, is visible. Fig. 2 illustrates that effect for three values of 9 which have 
different sign. The asymptotic, Markovian limit is achieved first for the subdiffusive case 6 — 0.4. 

For the case of the power law kernel we calculate the fractional diffusion coefficient by means of the Eq. ([5^ ; the 
quantity af'it) = F{t) is given by the Eq.dSO]). We obtain 



2ao7°r(l + a) 
TTfiT{a — 70; + 1) 



sin(/i7r/2)t 



2a— 7a — 1 



t 



[M(l-7)-9]/(M+e) 



(33) 



The diffusion properties of the system follow directly from the above formula. The influence of the parameter 9, which 
quantifies the structure of the medium, is similar as in Markovian case (20j : the larger is 9 the weaker is diffusion. For 
9 < 0, there is clearly the superdiffusion. For the positive 9, the diffusion becomes weaker with 9 and finally it turns 
to the subdiffusion; the critical value, which corresponds to the normal diffusion, is 9cr = /i(l —7). On the other hand, 
a < 9 < fi, there is a critical value of 7 which separates the superdiffusion from the subdiffusion: ^cr = 1 ^ ^//^- For 
9 > fi the motion is always subdiffusive. The parameter 7 measures the degree of the time nonlocality; it is the largest 
if 7 approaches 0. The diffusion speed grows if 7 diminishes because the system behaviour at large times becomes 
sensitive to the initial stages of the evolution when the distribution spreads rapidly. The latter conclusion shows 
that the memory can influence the diffusion in many ways: the non-Mar kovian CTRW predicts the weakening of the 
diffusion and it is just a consequence of the memory in the system [27| . However, in that case the time nonlocality 
invokes a trapping mechanism. 

Note that the above properties, in particular the presence of a transition from the subdiffusion to the superdiffusion 
when changing the parameters of the system, still hold if one considers some other fractional moment of order 5 < fi, 
instead of the renormalized moment A4^. 

For any kernel K^, except for the delta function and for the exponential kernel, the time evolution of the moment 
A4f^ is governed by the nonlocal equation (|30p and the diffusion properties follow from its solution. In fact, looking for 
the full solution may be avoided in some cases: the kind of diffusion is already determined by the sign of the function 
^{t) in the limit of long time. 



V. SUMMARY AND DISCUSSION 



We have studied the diffusion process for the non-Markovian systems with the position-dependent diffusion co- 
efficient, which involves Levy flights and then the variance of the corresponding probability distribution is infinite. 
The integral equation for that problem contains the fractional Riesz-Weyl operator and the time-dependent memory 
kernel; the diffusion coefficient depends on the position in the algebraic, scaling way. The equation has been solved 
in terms of the Fox functions in the limit of small wave numbers. We have demonstrated that this solution represents 
the Levy process for any memory kernel. The formal solution has been obtained in the closed form which involves 
the Laplace transform. The inversion of that transform may be a difficult task for the most of the kernels and then 
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numerical methods have to be apphed. Two forms of the kernel have been discussed in details: the exponential kernel, 
for which the problem resolves itself to the generalized telegrapher's equation, and power-law one which is equivalent 
to the fractional equation with both the Riesz-Weyl operator and the Riemann-Liouville fractional operator. For the 
exponential kernel, the memory initially slows down the spread of the distribution but asymptotically the solution 
converges to that of the Markovian equation. The case with the power-law kernel reveals much more interesting be- 
haviour. There is an interplay among all ingredients of the dynamics, in particular between the range of the memory 
7 and the inhomogeneity parameter 9, which can result in all kinds of the diffusion, both normal and anomalous. In 
order to make that classification possible, we have introduced the fractional diffusion coefficient, defined in terms of 
the renormalized moment of order /n. This coefficient allows us to maintain the standard terminology of the anomalous 
diffusion also for the Levy flights. 
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